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The photoproduction of beauty quarks in ep collisions has been measured using a data 
sample of 170 pb~ collected with the HI detector at HERA-II in the years 2006 and 
2007. Events with two jets and a muon in the final state were investigated, and beauty 
events were identified using the muon's relative transverse momentum to a jet and its 
impact parameter. Visible cross sections were measured differentially in the transverse 
momenta of the highest energy jet (p 1 -^ 1 ) and the muon (p^), the pseudorapidity of the 
muon (rf*) and of the photon's momentum fraction x° s entering the hard interaction. 
The measurements are found to be well described by QCD calculations at NLO. 

1 Introduction 

The production of beauty quarks in ep collisions has been investigated in considerable detail 
at HERA-I. In several analyses it was observed that measured cross sections where signif- 
icantly above the predictions from perturbative QCD calculations in next-to-leading order 
(NLO). The data from HERA-II with its larger statistics makes it possible to repeat these 
measurements with increased accuracy. 

This measurement Q] [2] follows closely a measurement performed with data from HERA- 
I [3], where beauty photoproduction events were investigated with two jets and a muon in the 
final state. The result of the HERA-I measurement was that NLO calculations describe the 
data reasonably well, except for the lowest bin of the muon and jet transverse momentum, 

and p!p , where the data were significantly above the predictions. Similar measurement 
have also been made by the ZEUS collaboration [5], in a slightly different phase space. Here, 
good agreement was found with QCD predictions, also at low jet and muon px- 

The measurement presented here uses the HERA-II data set to measure the same cross 
sections in the same visible range as in the previous HI publication with increased statistics 
and correspondingly smaller errors. The data was collected with the HI detector in the 
years 2006 and 2007, when HERA collided electron and positron beams of an energy of 
E c = 27.55 GcV with protons of E p — 920 GcV, and corresponds to an integrated luminosity 
of 170 pb" 1 . 



2 QCD Models 

The Monte Carlo generators PYTHIA 6.2 [5] and Cascade 2.0 [5] were used for the simulation 
of signal and background distributions. 

The PYTHIA event samples were generated with massless matrix elements, using the 
CTEQ6L [7] and SAS-1D [5] parton density sets for the proton and the photon, respectively. 
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The fragmentation of heavy quarks to hadrons was simulated using the Peterson fragmen- 
tation function with a parameter eb( c ) = 0.0069(0.058) for beauty (charm) quarks. More 
details on the parameter settings can be found in our previous publication [3]. 

For cross checks, additional Monte Carlo samples were generated using the Cascade 
program, which is based on fcx factorization and the CCFM evolution rather than collinear 
factorization and DGLAP evolution. The proton parton density set AO [9] is used for the 
unintcgratcd gluon density in the proton. 

Both programs, PYTHIA and Cascade, use Leading Order (LO), i.e. 0(a s ), QCD matrix 
elements for the hard scattering, augmented by parton showers that approximate the effect 
of additional multiple gluon emission. 

To compare the data to Next-to-Leading Order (NLO) QCD, i.e. 0(af), the FMNR 
program [TU] has been used, which is based on the NLO calculation by Nason, Dawson, and 
Ellis [II]. Here, the CTEQ5F4 [12] (GRV-G HO [13]) parton density set for the proton 
(photon) were employed. Hadronisation corrections were calculated using the PYTHIA 
program. 

The theory uncertainties were determined by varying the input beauty mass mh from the 
nominal value — 4.75 GcV up and down by 0.25 GeV. In addition, the renormalization 
and factorization scales /i r and //f were varied independently in the range /io/2 < /i r ,Mf ^ 
2 /^o, with the constraint 1/2 < /i r /jUf < 2, and no set to /io = v^Px + mi. For each bin on 
hadron level, the deviations due to the beauty mass variation and the largest deviation due 
to the scale variation in the upward and downward direction were determined and added in 
quadrature for the total model uncertainty, following the prescription in [14| p. 406]. 

3 Analysis Method 

Events were selected with the following experimental cuts: No electron candidate with an 
energy above E > 6 GcV is allowed to be found in the detector. The inelasticity y, measured 
from the hadronic activity in the detector, must lie in the range 0.2 < j/h < 0.8. Two jets 
with transverse momenta p^ 1 ^ > 7(6) GeV for the highest (second highest) px jet have 
to be found within the pseudorapitity range —2.5 < rf ct < 2.5. Jets were identified by the 
inclusive fc t jet algorithm [15] in the px recombination scheme, with a distance parameter 
R = 1.0. The jet algorithm was applied in the laboratory frame. A muon, identified in 
the central muon system, with transverse momentum pt^, > 2.5 GeV has to be found in a 
pseudorapidity range —0.55 < ?7 M < 1.1. The muon has to be associated to one of the two 
highest energy jets. 

To ensure a good muon reconstruction and suppress events with cosmic muons, additional 
quality cuts were applied to the muon track; in particular, the track must be associated to 
at least one hit in the central silicon tracker to ensure a precise measurement of the impact 
parameter 5. 

4 Cross Section Definition and Measurement 

We have measured the cross section of the process ep — » ebbX — » ejjfxX', i.e. beauty 
production with the formation of two jets and the subsequent decay of a beauty hadron to 
a muon. The muon may be produced by a direct semileptonic decay of a beauty hadron, 
from a cascade decay, where the charm hadron decays scmilcptonically, or from a J/ip or ip' 
decay. The muon must be associated to cither of the two highest px jets. 
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Figure 1 : Distributions of a) the impact parameter S of the muon track and b) the transverse 
muon momentum p™ 1 relative to the axis of the associated jet. Included in the figure are the 
estimated contributions of events arising from beauty quarks (dark grey line), charm quarks 
(black line) and light quarks (dotted line). 

The visible range is defined by Q 2 < lGeV 2 , 0.2 < y < 0.8, pt^ > 2.5 GeV, -0.55 < 
t/ 1 < 1.1 p J ^ tl(2) > 7(6) GeV, and -2.5 < rf ct1 ^ < 2.5. 

To extract the beauty fraction, two quantities were used: p™, which is the transverse 
momentum of the muon with respect to the axis of the most energetic jet, reconstructed 
without the muon four vector, and 6, which is the impact parameter of the muon with 
respect to the primary vertex of the event. The sign of 6 is determined in relation to the 
jet axis and defined such that it is positive for muons originating from a secondary vertex 
displaced along the jet direction. 

For each bin, the two-dimensional distribution of these quantities were fitted with three 
template distributions derived from Monte Carlo simulations, taking into account the statis- 
tical uncertainties from the data sample and the Monte Carlo templates [16] ; the templates 
were generated separately for events containing only light (u, d, and s) quarks, charm quarks, 
and beauty quarks, respectively. 

The result of the fit is the relative amount of beauty induced events /t, in each analysis 
bin, from which the observed number of beauty events N\> in the bin is calculated as iVb = 
fb ■ iVbinj where TVbm ^ s * ne total number of data events observed in the respective bin; the 
data were then corrected for effects of detector resolution by a matrix unfolding procedure, 
with a migration matrix determined from Monte Carlo simulations. 

Control Distributions 

Fig. [1] shows the distributions of the quantities used to extract the beauty fraction from 
the data: the impact parameter 6 of the muon track and the transverse muon momentum 
p™ 1 relative to the axis of the associated jet. Both quantities are described quite well by 
the Monte Carlo simulation. In particular, the impact parameter distribution, which is 
very sensitive to the detector resolution, is very well described, in the region 5 < that is 
dominated by resolution effects as well as in the region <5 > 0, which shows the tails due to 
long-lived particles from charm and beauty decays. 
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Systematic Uncertainties 

A number of sources of systematic errors were considered. The overall normalization un- 
certainty comprises uncertainties on the trigger efficiency, muon identification and track 
finding efficiencies, and the luminosity. Additional uncertainties that affect the data differ- 
ently in various bins arc: The impact parameter resolution, the reconstruction of the jet 
axis, the energy scale for hadrons of the calorimeter, the model uncertainties (estimated by 
using CASCADE [6] instead of PYTHIA 6.2 [5]), the uncertainty from the fragmentation 
process (estimated by using the Lund instead of the Peterson fragmentation function), the 
uncertainty from the fragmentation fractions of c and b quarks into hadrons, their branching 
ratios and lifetimes, and the uncertainty on the modelling of ir and K inflight decays. The 
resulting systematic uncertainty is 12%. 

5 Results 

The total cross section for the process ep — > ebbX — > ejjfj,X' in the visible range given above 
has been measured to be 

CTvis (ep -> ebbX -> ejjfxX') = 31.4 ± 1.3(stat.) ± 3.8(sys<.) pb. 

This result is somewhat lower than the published result [3j from HERA-I, which is 
(7 V i S = 38.4 ± 3.4(siai.) ± 5A(syst.) pb, but compatible within errors. In comparison, the 
FMNR calculation yields <r V is = 25. pb, in agreement with the data, and the PYTHIA 
prediction is a v - ls = 21.7 pb. 

Differential cross sections were also measured as a function of the following quantities: 
The transverse momentum of the muon p T , the transverse momentum pi^ 1 of the highest p^ 
jet, the pseudorapidity rf 1 of the muon, x° hs , the momentum fraction of the photon entering 
the hard interaction, and <5</>j C ts, the difference in azimuthal angle between the two jets. 

The differential cross sections also tend to be lower than the HERA-I results, as shown 
in Fig. [21 The largest discrepancies arc observed for the differential measurements in the 
lowest bins of p T and p 3 ^ 1 . For these bins, the ratio between the measurements is about 
2.5 a below unity if systematic uncertainties that do not cancel between both measurements 
are taken into account. It has been checked that this discrepancy is not caused by differences 
in the analysis method between the HERA-I and HERA-II analyses; we therefore attribute 
the difference to a statistical fluctuation. 

Fig. [2k) shows the difference in azimuthal angle, <5</>j e ts, between the two jets. In leading 
order QCD the two outgoing quarks must be exactly opposite in azimuthal angle, correspond- 
ing to Scj)j e ts = 180°. Thus, values of 50j o ts substantially lower than 180° are indicative of 
the presence of further final state gluons, and therefore this quantity is sensitive to the de- 
scription of gluon emission. PYTHIA, which employs parton showers to simulate the effect 
of multiple gluon emission, describes the shape of this observable reasonably well, as does 
the fixed order calculation of FMNR, which allows at most one hard gluon in the final state. 

Overall, the data are reasonably well described in shape by the predictions from PYTHIA, 
but lie approximatly a factor 1.4 above the PYTHIA prediction. The NLO predictions, 
derived with the FMNR program, also lie systematically below the data, but also describe 
the differential distributions well in shape. In particular, the deficiency in the lowest bins of 
p T and p!^ is not substantially larger than in the other bins, in contrast to the findings of 
the HERA-I analysis. This observation agrees with the results from ZEUS [4j. 
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Figure 2: Differential cross sections for the photoproduction process ep — > ebbX — ► ejjfiX in 
the kinematic range Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p£ > 2.5 GeV, 0.55 < jf < 1.1, p J ^ tl(2) > 
7(6) GeV and |yyj ctl ( 2 )| < 2.5. The cross sections are shown as functions of a) the muon 
pseudorapidity r]^, b) the muon transverse momentum p!^, c) the jet transverse momentum 
p^* 1 °f the highest transverse momentum jet, d) the photon's momentum fraction i° bs 
entering the hard interaction, and e) the azimuthal angle difference S<j)j e t B between the jets. 
The inner error bars show the statistical error, the outer error bars represent the statistical 
and systematic uncertainties added in quadrature. The NLO QCD predictions are corrected 
to the hadron level (solid line) using the PYTHIA generator. The shaded band around 
the hadron level prediction indicates the systematic uncertainties as estimated from scale 
variations (see text). Predictions from the Monte Carlo generator PYTHIA (dotted line) 
are also shown. 
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6 Conclusions 



We have performed a measurement of the photoproduction of beauty quarks, using events 
where at least one beauty hadron decays with a muon in the final state, and two jets are 
visible in the detector, in the phase space region defined by Q 2 < 1 GeV 2 , 0.2 < y < 0.8, p^ > 
2.5 GeV, -0.55 < if < l.l,p^ tl{2) > 7(6) GeV, and -2.5 < if ct1 ^ < 2.5. 

The visible cross section has been measured to be er v i s (ep — > ebbX — > ejj[iX') = 31.4 ± 
1.3(stat.) ± 3.8(syst.) pb. A NLO QCD calculation is in agreement with this measurement 
within the theoretical uncertainties. 

Differential cross sections have been measured as function of the observables ?7 M , pt^, 
p!^ , x° bs , and S4>j ets . The shape of these distributions is reasonably well described by the 
NLO QCD calculation as well as the PYTHIA LO Monte Carlo program. 

At low values of pi^ and p 3 ^ 1 , the new measurement lies lower than the previous HERA-l 
measurement published by HI, and is thus better described by the NLO predictions than 
the previous measurement. 
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